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Abstract 

Using 58xl0 6 J/ijj events collected with the Beijing Spectrometer (BESII) at the Beijing 
Electron Positron Collider (BEPC), the branching fractions of J ftp to ppr] and pprj' are de- 
termined. The ratio ^(jf-ip^pp) obtained by this analysis agrees with expectations based on 
soft-pion theorem calculations. 

PACS numbers: 
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I. INTRODUCTION 



The J/ip meson has hadronic, electromagnetic, and radiative decays to light hadrons, 
and a radiative transition to the i] c . In Ref. [1], direct hadronic, electromagnetic and 
radiative decays are estimated to account for 69.2%, 13.4%, and 4.3%, respectively, of 
all J/ip decays. However, individual exclusive J /if) decays are more difficult to analyze 
quantitatively in QCD. To date, two-body decay modes such as J ftp — > B 8 B$ or PgVg, 
where B s , P 9 and V 9 refer to baryon octet, pseudoscalar nonet, and vector nonet particle, 
respectively, have been studied with some success using an effective model, and other 
similar methods [2]. 

Studies of three-body decays of J /if) are a natural extension of studies of two-body 
decays. Since most J / if) decays proceed via two-body intermediate states, including wide 
resonances, it is hard to experimentally extract the non-resonant three-body contribu- 
tion [3]. Specific models based on proton and N* pole diagrams have been introduced to 
deal with these problems [2] . In the calculation, the soft-pion theorem [4] has been applied 
to the decay J /if) — > ppn successfully. This method has also been used for J /if) — > pprj 
and J /if) — > pprj' decays [2] . 

This paper reports measurement of the branching fractions for pprj and pprj' , and tests 
of the soft-pion theorem for J /if) — > pprj, which states [2]: 

Y{J/if) — > pprj) 



r(j/tp^pp) 



~ 0.64 ±0.52. 



II. THE BES DETECTOR AND MONTE CARLO SIMULATION 



BESII is a conventional solenoidal magnet detector that is described in detail in 
Refs. [5]. A 12-layer vertex chamber (VC) surrounding the beam pipe provides trigger 
and track information. A forty-layer main drift chamber (MDC), located radially outside 
the VC, provides trajectory and energy loss (dE/dx) information for tracks over 85% of 
the total solid angle. The momentum resolution is a v /p = 0.017-^/1 + p 2 (p in GeV/c), 
and the dE/dx resolution for hadron tracks is ~ 8%. An array of 48 scintillation counters 
surrounding the MDC measures the time-of-flight (TOF) of tracks with a resolution of 
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~ 200 ps for hadrons. Radially outside the TOF system is a 12 radiation length, lead-gas 
barrel shower counter (BSC). This measures the energies of electrons and photons over 
~ 80% of the total solid angle with an energy resolution of a E /E = 22%/ ^fE (E in GeV). 
Outside of the solenoidal coil, which provides a 0.4 Tesla magnetic field over the tracking 
volume, is an iron flux return that is instrumented with three double layers of counters 
that identify muons of momentum greater than 0.5 GeV/c. 

In the analysis, a GEANT3-based Monte Carlo (MC) simulation program (SIMBES) [6] 
with detailed consideration of detector performance is used. The consistency between data 
and MC has been validated using many high purity physics channels [7]. 

In this analysis, the detection efficiency for each decay mode is determined by a MC 
simulation that takes into account the angular distributions. For J/ip — > ppi], the angle 
(9) between the directions of e + and p in the laboratory frame is generated according to 
1 + a ■ cos 2 6 distribution, where a is obtained by fitting the data from J/ip — > pprj. A 
uniform phase space distribution is used for J ftp decaying into pprj' . 

III. GENERAL SELECTION CRITERIA 

Candidate events are required to satisfy the following common selection criteria: 

A. Charged track selection 

Each charged track must: (1) have a good helix fit in order to ensure a correct error 
matrix in the kinematic fit; (2) originate from the interaction region, ^V 2 + V 2 < 2 cm 
and \V Z \ < 20 cm, where V x , V y , and V z are the x, y and z coordinates of the point of 
closest approach of the track to the beam axis; (3) have a transverse momentum greater 
than 70 MeV/c; and (4) have | cos#| < 0.80, where 6 is the polar angle of the track. 
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B. Photon selection 

A neutral cluster in the BSC is assumed to be a photon candidate if the following 
requirements are satisfied: (1) the energy deposited in the BSC is greater than 0.05 GeV; 
(2) energy is deposited in more than one layer; (3) the angle between the direction of 
photon emission and the direction of shower development is less than 30°; and (4) the 
angle between the photon and the nearest charged track is greater than 5° (if the charged 
track is an antiproton, the angle is required to be great than 25°). 

C. Particle Identification (PID) 

For each charged track in an event, Xpid(^) * s determined using both dE/dx and TOF 
information: 

X 2 pi D (i) = X^E/dx^) + Xtof(i), 

where % corresponds to the particle hypothesis. A charged track is identified as a pion if 
Xpid f° r the 7r hypothesis is less than those for the kaon and proton hypotheses. For p or 
p identification, the same method is used. In this analysis, all charged tracks are required 
to be positively identified. 

IV. ANALYSIS OF J/rj) ppr) 

The decay modes for the J/ijj — > ppr] measurement are f] — > 77 and 77 — > 7r + 7r~7r°. 
The use of different decay modes allows us to cross check our measurements, as well as 
to obtain higher statistical precision. 

A. r\ —> 77 

Events with two charged tracks and two photons are selected. A four-constraint (4C) 
kinematic fit is performed to the hypothesis J/tp — > M077. For events with more than two 
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FIG. 1: The two-photon invariant mass distribution for J/ip — > ppj^y candidate events. Data 
are represented by rectangles; the error bars are too small to be seen. The curves are the results 
of the fit described in the text. The shaded part is background from MC simulation. 



photons, all combinations are tried, and the combination with the smallest x 2 is retained. 
X^-ypp is required to be less than 20. 

The 77 invariant mass (m 77 ) distribution for selected events is shown in Fig. 1. A peak 
around the r\ mass is evident. The curves in the figure indicate the best fit to the signal 
and background. The shaded part is the background estimated from a MC simulation of 
inclusive J/ip events [8]. The main background comes from J/ip — > ppir°-ir and By 
fitting the r\ signal with a MC-simulated signal histogram plus a third order polynomial 
background function, the number of r\ signal events is determined to be (12220 ± 149). 

For the signal MC simulation, the events are generated with a proton angle distribution 
of 1 + a cos 2 9, where a is taken to be -0.6185 in order to describe the data. In the decay, 
intermediate resonances, N(1440), N(1535), N(1650), and N(1800) and antiparticles, with 
fractional contribution of (8 ±4)%, (56 ±15)%, (24^ 5 )%, and (12 ±7)% [9], respectively, 
are included. The resulting detection efficiency for J/ip — > ppi] (rj — > 77) is determined 
to be 28.70%. The ppi] branching fraction, calculated using 

B(J/ip -> ppi]) = — — , 

e • N J/lP ■ B(r) -> 2 7 ) • /i 

is (1.93 ±0.02) x 10~ 3 , where the error is statistical only. Here N obs represents the number 
of observed events, e is the detection efficiency for J/ip — > ppr\{j] — > 27), f\ = 0.9739 is the 
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efficiency correction factor (see Section VI), and Nj/^ is the total number of J/ip events. 



400 




FIG. 2: The 7r + 7r~7r invariant mass distribution for J/ip — > pptt + tt~tt° candidate events. The 
curves are results of the fit described in the text. The shaded part is background from MC 
simulation. 



B. 7) —* 7T+7T 7T° 

Similar to the above analysis, events with four charged tracks and two photons are 
selected. A 4C kinematic fit is performed to the J/ip — > ppix+ix"^ hypothesis, and 
the X^ppn+n- value is required to be less than 20. In order to suppress multi-photon 
backgrounds, the number of photons is required to be two. The invariant mass of the 77 
is required to be between 0.095 and 0.175 GeV/c 2 . 

The 7r + 7r~7r° invariant mass (m^+^-^o) distribution is shown in Fig. 2, where a peak at 
the 7] mass is observed. The curves in the figure are the results of a fit to the signal and 
background. The shaded part is background estimated from MC simulation of inclusive 
J/ip decay events [8]. Here the main background comes from J/xjj — > ppn + 7T~7T° and 
pp7r + 7r~7 decays. By fitting the distribution with a MC simulated signal histogram plus 
a third order polynomial background function, (954 ± 45) signal events are obtained. 
Similar to the rj — > 27 decay, contributions from the baryon excited states N(1440), 
N(1535), N(1650), and N(1800), as well as their anti-particles [9], are considered. The 
detection efficiency of J ftp — > pprj (rj — > -k + it~it , tt° — > 77) is determined to be 4.20%. 
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The branching fraction is determined from the calculation 

B(J/ip^ppr]) = — _ 

e • Nj/^ ■ B(T] -> 7T + 7T 7T°) 
1 

5(vr ^77)-/ 2 ' 

where / 2 = 0.9582 is a correction factor for the efficiency that is described below in 
Section VI. We determine a branching fraction for J/ijj — > ppr/ of (1.83 ± 0.09) x 10~ 3 , 
where the error is statistical only. 
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FIG. 3: The 7r + 7r _ ?7 invariant mass distribution for J/ip — > ppir + ir~ri candidate events. The 
curves are results of the fit described in the text. The shaded part is background from MC 
simulation. 



V. ANALYSIS OF J/V> -> pprj' 

There are three main decay modes of the ?/: rf — > 7t + tt~i], r)' — > 7p° and 77' — > 7r 7r°77. 
Here the first two decay modes are used. 



A. 77' — ► 7r+7r 77, 77 — ^ 77 

In the search for 7/ — > n + n^r] decays, events with four charged tracks and two photons 
are selected. A five-constraint (5C) kinematic fit is performed to the hypothesis of J/ip — > 
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ppTT + 7T 77, in which the 27 invariant mass is constrained to equal the r) mass, and the 
y 2 _ + _ value is required to be less than 20. 

77PP 71 " 

The 7r + iT~ri invariant mass (m 7r + 7r - r? ) distribution for events that survive the selection 
criteria is shown in Fig. 3. A clear rf signal is observed. The curves in the figure are 
the best fit to the signal and background. The shaded part is background estimated 
from MC simulation of inclusive J /if) decay events [8]. The main background comes from 
J/ip — > A + A~r], and A°A°r] decays. By fitting the distribution with a MC simulated 
signal histogram plus a third order polynomial background function, a signal yield of 
(65 ± 12) events is observed. According to a MC simulation, in which the events are 
generated with uniform phase space, the detection efficiency of J/tp — > ppi]' {rf — > n + 7Y^r], 
V ~^ 77) is 3.38%. The effect of intermediate resonances is considered as a source of 
systematic error. Using 

B{J/i)^ppq') = — — - 

e • Nj/^ ■ B{rf -> tt+tt rj) 

1 

B(n -> 77) • / 3 

with / 3 = 0.8228 being the efficiency correction factor (see Section VI), we determine the 
branching fraction for J/ip — > ppv( to be (2.31 ±0.43) x 10 -4 , where the error is statistical 
only. 
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FIG. 4: The 77r + 7r~ invariant mass distribution for J/ip — > ppjir + iT~ candidate events. The 
curves are results of the fit described in the text. The shaded part is background from MC 
simulation. 
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B. T)' — ► 7p°, p° — > 7T+7T- 

In order to select p' — > 7p°, a 4C kinematic fit is performed under the hypothesis of 
J/ip —> pp7r + 7r~7. The x^ppTT+v,-- value is required to be less than 20. To ensure the events 
are from 7p°, a \m n + n - — m p \ < 0.20 GeV/c 2 requirement is imposed, where m n + w - is the 
7r + 7r _ invariant mass, and m p is the p mass. In order to exclude the background from 
J/ip — > pp7i + 7i~ , it is required that the invariant mass of the four charged tracks is less 
than 3.02 GeV/c 2 . 

The 7p° invariant mass (m 7p o) distribution for selected events, where a clear p' signal is 
observed, is shown in Fig. 4. The curves in the figure are the best fit to the signal and back- 
ground. The shaded part is the background estimated from MC simulation of inclusive 
J/ip decay events [8]. The main background comes from J/t/j — > pp7r + 7r~7, A ++ A it 
and pp7T + 7r~7r° decays. By fitting the m 17r + % - distribution with a MC simulated signal 
shape and a third order polynomial background function, we determine the number of r/' 
signal events to be (200 ±29). The detection efficiency for J/ip ppr]' (7/ — > 7p°) is deter- 
mined to be 7.48%, assuming phase space production, where the 7r + 7r~ mass distribution 
is generated according to measurements fromJ/ip — > (f)r)',r)' — > •yn + ii~ [10]. Using 

B(J/ip^ppr)') = — 

e • iVj /v , • B(rf -> 7p°) • / 4 

with the /4 correction factor of 0.8522 (see Section VI). The resulting branching fraction 
for J/ip — > ppp' is (1.85 ± 0.27) x 10~ 4 , where the error is statistical only. 

TABLE I: Numbers used in the calculations of branching fractions and results. 

Decay mode N b s e(%) Branching fraction 

J/iP -► pprj, 7] -> 77 12220 ± 149 28.70 B{ J/%b -► ppr/) = (1.93 ± 0.02 ± 0.18) x 10" 3 
J/ip -»• ppr/, r? tt+tt-tt 954 ± 45 4.20 £( J/^ ppq) = (1.83 ± 0.09 ± 0.24) x 10~ 3 
J/tp -► ppr/', r/ -► vr+Tr-r/ 65 ± 12 3.38 £( J/Y> -► ppr/') = (2.31 ± 0.43 ± 0.34) x 10~ 4 
J/V> PPV, r/' -► 7/9° 200 ± 29 7.48 B(J/ip -► ppr/') = (1.85 ± 0.27 ± 0.31) x 10~ 4 
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VI. SYSTEMATIC ERRORS 



In our analysis, the systematic errors on the branching fractions come from the uncer- 
tainties in the MDC tracking, photon efficiency, particle identification, photon identifica- 
tion, kinematic fit, background shapes, hadronic interaction model, intermediate decay 
branching fraction, the ir° and p selection requirements, intermediate resonance states, 
and the total number of J/ip events. The errors from the different sources are listed in 
Table II. 

The MDC tracking efficiency has been measured using J ftp — > pir, AA, and ip(2S) — > 
n + ii^J/ip, J /if) to p + p~. The MC simulation agrees with data within 1 to 2% for each 
charged track [7]. Thus 4% is regarded as the systematic error for the two charged-track 
mode, and 8% for the four charged-track final states. 

The photon detection efficiency has been studied using a sample of J/ip — > pir [7] 
decays; the difference between data and MC simulation is about 2% for each photon. In 
this analysis, 2% is included in the systematic error for one-photon modes and 4% for 
two-photon modes. 

The charged pion PID efficiency has been studied using J/ip — > pir decays [7]. The PID 
efficiency from data is in good agreement with that from MC simulation with an average 
difference that is less than 1% for each charged pion. Here 2% is taken as the systematic 
error for identifying two pions. 

The proton PID efficiencies have been studied using J/ip — > ppn + iT~ decays. The main 
difference between data and MC simulation occurs for tracks with momentum less than 
0.35 GeV/c. We determine a weighting factor for identifying a proton or anti-proton as a 
function of momentum from studies of the J/ip — > ppTi + n~ channel. After considering the 
weight of each particle in an event, the difference between data and MC is determined to 
be g£ = 0.9739±0.0078 for 77 -> 2 7 , 0.9582±0.0199 for rj -> 3tt, 0.8228±0.0211 for rj' -> 
n+n-r], and 0.8522 ± 0.0140 for rf -> 7p°. We take f x = 0.9739, f 2 = 0.9582, f 3 = 0.8228, 
and fi = 0.8522 as efficiency correction factors for the corresponding decay channel, and 
0.8%, 2.1%, 2.6%, and 1.6% are taken as the errors associated with identifying protons 
and anti-protons, respectively. The PID systematic errors for the four decay modes are 
listed in Table II. 
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For the systematic error of photon ID, which arises mainly from the simulation of fake 
photons, pp and J/ip — > ppn + n^ data samples were selected and 10 5 simulated pp and 
J /if) — > ppix + ir~ events were generated, with real and fake photons. The decay mode 
J/ip ^ pp is used for the photon ID systematic error of J/ip — > pprj (rj — > 27), and 
the decay mode J/^ — > pp7r + 7r~ for J/^ — > pp7/ (77 — > 37r) and J/^ — > pp7/. From 
the decay mode J/^ — > pp, the fake photon differences between data and MC is about 
2.0%, while for the decay mode J/V — > ppn + ir~, the difference is 1.6%. Here 2.0% is 
taken as the systematic error associated with photon ID for the decay mode determined 
to be J /if) — > pp7/ (77 — > 77), and 1.6% for the decay modes J/^ — > pp7/ (77 — > 37r) and 
J/V> — > pp7/. 

In Ref. [11], the uncertainty of the 4C kinematic fit is 4%, which we include here in 
the systematic error. The uncertainty of the 5C kinematic fit is 4.1% in Ref. [12]. Here 
we conservatively take 5% as the systematic error from the 5C kinematic fit for the decay 
mode 7/ — > 7r + 7r~T7. 

The systematic errors of the background uncertainty is obtained by changing the range 
of the fit and varying the order of the polynomial background. The errors range from 
0.8% to 7.3% in all decay modes (see Table II for detail). 

There are two models, FLUKA and GCALOR, used for simulating hadronic interac- 
tions; the different models lead to different detection efficiencies. The difference between 
them is regarded as a systematic error. For the decay J ftp — > pprj (rj — > 27), the difference 
is very small and negligible. For the other decay modes, it is about 1.4% for J/ip — > pprj 
(rj — > 7T + 7r~7r°), J/xjj — > pprj' (rj' — > 7i + n^ri), and 5.2% for J/ip — > pprj' (rj' — > 7p°). 

The branching fractions for the decays 7r° — > 27, 77 — > 27, 77 — > 7r + 7r~7r°, 77' — > 7i + 7T~i], 
and 77' — > 7p are taken from the PDG [13]. The errors on these branching fractions are 
systematic errors in our measurements. 

For the 77 — > 3ir mode, the n° mass is required to satisfy |M 77 — M n o\ < 0.04 GeV/c 2 . 
To study the systematic error associated with this requirement, 7r° samples are selected 
and simulated using J/ip — > p7r, and the data and MC efficiencies in the 3cr signal region 
are compared with using the requirement or not, the difference is about 1%. Here it is 
taken as the systematic error caused by the n° requirement. 
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For the rj' — > 7p mode, we require that |M 7r + 7r - — M p | < 0.20 GeV/c 2 . According to 
Ref. [14], the uncertainty associated with this requirement is 5.9%. Here we take this as 
the systematic error for the p mass requirement. 

In the signal MC simulation, we assume the presence of N(1440), N(1535), N(1650), 
and N(1800) in the pprj channel. If some of these resonances are not included, the efficiency 
of this channel changes. These differences are taken as the systematic error associated 
with possible intermediate states. The total systematic error associated with this is taken 
as the sum added in quadrature. For the decay modes with an 7/, we take the difference in 
efficiency determined assuming the decay proceeds via an intermediate iV(2090) resonance 
compared with phase space generation as the systematic error (see Table II for detail). 

The uncertainty of the total number of J/ip events is 4.7% [15]. Combining all errors 
in quadrature gives total systematic errors of 9.3% for rj — > 77, 12.9% for rj — > it + it~it , 
14.8% for 7/ — > 7r + 7r~?7, and 16.6% for 7/ — > 7p. 



VII. RESULTS 

Table I shows the branching fractions of the two channels into their different decay 
modes; the first error is statistical and the second is systematic. The results for the 
different decay modes in the same channel are consistent within errors and are combined 
after taking out the common systematic errors (8.37 % for the rj mode and 10.8% for the 
77' mode): 

Br(J/i> -> pprj) = (1.91 ± 0.17) x 10~ 3 , 
Br(J/iP -> pprj') = (2.00 ± 0.36) x 10~ 4 . 

In comparison with previous measurements of J/ip — > pprj and J/ip — > pprj' , the present 
results are of higher precision. 

Using the result of J/ip — > pprj from this analysis and that of J/tp — > pp in Ref. [16], 
we determine: 

m^M}! ^ q.85 ±0 .08, 

r( j/V - pp) 
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TABLE II: Summary of systematic errors; "-" means no contribution. 



Sources 




Relative error (%) 






Decay modes 


rj -> 27 


77 — > 7T + 7T~7r 7/' - 


-> 7T+7T 


'rj rj' — > 7/j° 


MDC tracking 


4 


8 


8 


8 


Photon detection efficiency 


4 


4 


4 


2 


Particle ID 


~1 


4.1 


4.6 


3.6 


Photon ID 


2.0 


1.6 


1.6 


1.6 


Kinematic fit 


4.0 


4.0 


5.0 


4.0 


Background uncertainty 


~1 


3.1 


7.3 


5.8 


Hadronic Interaction Model 


~0 


1.4 


1.4 


5.2 


Xlltcl IIlcU-lcxLc LLcCdj J_)l . 171. 


~ 1 

^ 1 


1 9 

1 . Zj 


3 1 

O. -L 


3 1 

O. -L 


7T° selection 




~1 






p selection 








5.9 


Intermediate resonances 


3.0 


4.0 


2.0 


7.1 


Number of J/ip events 


4.7 


4.7 


4.7 


4.7 


Total systematic error 


9.3 


12.9 


14.8 


16.6 



This is consistent with the calculation based on the soft-pion theorem, and indicates that 
the contribution of iV*- pole diagram is dominant for the J/i/j — > ppi] mode. 
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